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Abstract 


At  the  present  time  RANS  calculations  are  widely  used  to  determine  the  time  averaged  flow  over  bodies 
of  arbitrary  shape.  However  they  do  not  provide  the  time  dependent  or  spectral  information  required  for 
the  analysis  of  ship  hull  vibration,  radiated  noise  and  sonar  performance.  The  objective  of  this  study  is  to 
calculate  space  time  velocity  correlation  functions  and  surface  pressure  wavenumber  spectra  in  turbulent 
shear  flows  using  information  from  time  averaged  RANS  calculations,  so  that  the  details  of  high  fidelity 
time  averaged  flow  can  be  applied  to  hydroacoustic  and  structural  vibration  problems. 

Recent  advances  in  numerical  techniques  have  resulted  in  efficient  methods  for  calculating  the  Reynolds 
average  statistics  of  the  flow  around  most  bodies.  These  calculations  provide  the  mean,  time  invariant 
flow  as  well  as  the  local  turbulent  kinetic  energy  and  dissipation.  However  they  provide  no  information 
on  the  higher  order  statistics  or  spectral  content  of  the  unsteady  part  of  the  flow.  For  many  applications, 
especially  those  involving  sound  radiation,  the  local  turbulent  kinetic  energy  docs  not  provide  the 
information  that  is  required.  Ideally,  a  procedure  for  estimating  the  two  point  correlation  function  of  the 
velocity  fluctuations  throughout  the  flow  is  needed. 

In  a  recent  study  (Glegg  et  tf/(2010)),  a  method  for  calculating  the  surface  pressure  spectra  below  a 
turbulent  boundary  layer  was  developed  which  is  based  on  the  solution  to  the  unsteady  Euler  equations 
and  the  Reynolds  averaged  Navier  Stokes  equations.  It  was  shown  that  the  unsteady  velocity  in  a  two 
dimensional  linearized  shear  flow  could  be  modeled  using  a  distribution  of  vortex  sheets.  Each  vortex 
sheet  is  convected  at  the  local  mean  flow  speed  and  its  distribution  can  be  used  to  specify  all  three 
components  of  the  turbulent  flow.  The  unsteady  pressure,  the  two  point  velocity  correlation  functions,  the 
turbulent  kinetic  energy  and  surface  pressure  wavenumber  spectra  was  obtained  directly  from  this 
equation  in  terms  of  the  mean  square  value  of  the  vortex  sheet  strength.  Since  RANS  calculations  give  the 
distribution  of  turbulent  kinetic  energy  in  the  flow,  it  was  possible  to  invert  the  spatial  distribution  of 
turbulent  kinetic  energy  to  obtain  the  vortex  sheet  strength,  and  use  the  result  to  calculate  all  the  statistical 
details  of  the  unsteady  flow.  However,  in  this  approach  the  specification  of  a  turbulent  energy  spectrum 
that  scales  on  the  dissipation  length  scale  was  required  and  this  is  difficult  to  model  universally. 

In  the  first  part  of  this  study  (see  Part  la)  the  linearized  approach  described  above  was  extended  to 
turbulent  flows  in  two  dimensional  wakes  described  by  the  Orr  Sommerfeld  equations.  Wake  flows  are 
more  challenging  than  two  dimensional  boundary  layers  because  the  flow  is  continuously  evolving  in  the 
downstream  direction  in  a  self  similar  manner.  The  theory  was  modified  to  allow  for  self  similarity  and 
the  results  were  compared  to  experimental  data  of  the  two  point  velocity  correlation  functions  in  turbulent 
wakes  (Devenport  et  al,  2001)  with  some  success. 

This  approach  is  then  used  to  predict  rotor  stator  interaction  noise  that  is  caused  by  the  wakes  from 
upstream  blades  impinging  on  downstream  stator  vanes.  The  turbulence  modeling  developed  in  Part  la  is 
used  in  Part  lb  to  specify  the  turbulence  in  the  inflow  to  the  stators.  Acoustic  predictions  show  how  the 
sound  radiation  is  affected  by  the  cutting  angle  of  the  wake  relative  to  the  leading  edge  of  the  blades. 

In  the  second  part  of  this  report  (Part  II)  the  acoustic  radiation  from  a  turbulent  flow  over  a  step  is 
analyzed.  This  work  helps  to  explain  some  of  the  differences  between  the  measurements  by  Catlett  (2010) 
in  a  companion  study  at  Virginia  Tech  and  the  LES  based  predictions  presented  by  Ji  and  Wang(2010). 
The  previous  studies  of  the  problem  had  used  a  low  frequency  approximation  for  the  Greens  function. 

The  work  presented  here  shows  that,  at  the  frequencies  and  flow  Mach  numbers  considered  by 
Catlett(2010),  the  low  frequency  approximation  leads  to  incorrect  predictions  of  the  far  field  sound.  The 
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exact  Greens  function  for  a  step  is  given  in  Part  II  of  this  report  and  the  importance  of  the  diffracted  field 
is  identified. 

Part  III  of  this  report  focuses  on  flat  plate  boundary  layers  and  the  flow  downstream  of  a  forward  facing 
step  as  measured  Dr.  William  Devenport  in  a  companion  study  at  Virginia  Tech.  The  previous  analyses 
used  the  inviscid  form  of  the  linearized  equations  of  motion  to  model  the  turbulent  flow  and  in  the  last 
year  the  modeling  has  been  improved  to  include  viscous  terms.  In  addition  the  non  linear  interactions  are 
included  as  source  terms  for  fully  developed  turbulence.  The  results  have  been  compared  to  the 
experimental  measurements  reported  by  Awasthi  et  al  (201 1)  and  predictions  are  made  of  the  velocity 
spectra  and  the  surface  pressure  spectra. 


3 


Table  of  Contents 


Page 


Part  la:  Turbulence  Modeling  for  Plane  Wake  5 

Part  lb:  Turbulence  Modeling  for  Rotor  Stator  Interaction  Noise  10 

Part  2:  The  Tailored  Greens  Function  for  a  Step  2 1 

Part  3:  Predictions  of  Surface  Pressure  Spectra  for  a  Turbulent  Boundary  Layer  37 


4 


Part  la:Turbulence  Modeling  for  a  Plane  Wake 


1.1  Introduction 


In  a  recent  study  (Glegg  et  67/(20 1 0)),  a  method  for  calculating  the  surface  pressure  spectra  below  a 
turbulent  boundary  layer  was  developed  which  is  based  on  the  solution  to  the  unsteady  Euler  equations 
and  the  Reynolds  averaged  Navier  Stokes  equations.  It  was  shown  that  the  unsteady  velocity  in  a  two 
dimensional  linearized  shear  flow  could  be  modeled  using  a  distribution  of  vortex  sheets.  Each  vortex 
sheet  is  convected  at  the  local  mean  flow  speed  and  its  distribution  can  be  used  to  specify  all  three 
components  of  the  turbulent  flow.  The  unsteady  pressure,  the  two  point  velocity  correlation  functions,  the 
turbulent  kinetic  energy  and  surface  pressure  wavenumber  spectra  was  obtained  directly  from  this 
equation  in  terms  of  the  mean  square  value  of  the  vortex  sheet  strength.  Since  RANS  calculations  give  the 
distribution  of  turbulent  kinetic  energy  in  the  flow,  it  was  possible  to  invert  the  spatial  distribution  of 
turbulent  kinetic  energy  to  obtain  the  vortex  sheet  strength,  and  use  the  result  to  calculate  all  the  statistical 
details  of  the  unsteady  flow.  However,  in  this  approach  the  specification  of  a  turbulent  energy  spectrum 
that  scales  on  the  dissipation  length  scale  was  required  and  this  is  difficult  to  model  universally. 

In  the  first  part  of  this  study  the  approach  described  above  was  extended  to  turbulent  flows  in  tw  o 
dimensional  wakes.  Wake  flows  are  more  challenging  than  two  dimensional  boundary  layers  because  the 
flow  is  continuously  evolving  in  the  downstream  direction  in  a  self  similar  manner.  The  theoretical 
approach  was  modified  to  allow  for  self  similarity  and  the  results  w  ere  compared  to  experimental  data  of 
the  two  point  velocity  correlation  functions  in  turbulent  wakes  (Devenport  et  al,  2001)  with  some  success. 

1.2  Theory  for  the  turbulence  of  a  wake  flow 


Wc  will  consider  small  velocity  perturbations  (w,v,w)  in  a  parallel  shear  flow  defined,  in  eatesian 
coordinates  (x,y,z),  by  the  mean  velocity  U(y)  in  the  jc  direction  (see  figure  1).  The  linearized 
form  of  the  Navier  Stokes  equations  and  the  continuity  equation  can  be  used  to  obtain  the  Orr 
Sommerfcld  equation  for  the  upwash  component  of  the  velocity  (see  for  example  Hallbaek  1996) 


d'U  <&' 


-d-v4v  =  0 

R 

e 


(1) 


For  a  turbulent  flow  vvhieh  is  time  stationary,  but  may  not  be  homogeneous  in  the  jc  direction,  we 
will  consider  the  solution  of  this  equation  in  the  form 

(2) 


\ 1  (n)s  \ 

V  =  2rf  anV  (y)e 

/!=] 


J  I  vwOO  \2  dy  =  L(x) 


where  an  are  the  eigenvalues  of  the  solution  and  are  assumed  to  be  slowly  varying  in  the  flow 
direction,  so  the  flow  is  at  least  locally  homogeneous. 

The  other  components  of  the  velocity  can  be  obtained  by  considering  the  vortieity  in  the  y 
direction  given  by 
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(3) 


_  dit  dw 

T]~lk~!k 

Combining  this  with  the  continuity  equation  (Hallbaek,1996)  then  gives  the  mode  shapes  for  the 
other  two  velocity  components  as 

(4) 

«'” *>'" =4(P^+«,n'h') 

K  °y  K  dy 

where  A]  =a~,  +  [l1 .  To  bring  the  system  to  closure  we  use  the  vorticity  equation  for  the  r] 
component  which  gives 

(5) 

Dr] 1_V2  dU  a- 

D(  R  ^  dy  d: 


and  at  high  Reynolds  number  this  gives 

ayn'(y)e-“aita"'-‘*  =0 

ck  oy}  ,/=) 


(6) 


The  solution  to  this  equation  is  then  readily  obtained  by  assuming  that  the  vortiety  is  initially 
zero  and  gives  the  modes  of  the  vortieity  in  terms  of  the  upwash  velocity  modes 

(7) 

n=  1 

n(n)(y)  =  i 


Pv(n'(y)8n((0,x,y)  8n((0,x,y)=[— - 


U  dv 


i(an-(o/U)  x 


Hence  equations  (2),  (4)  and  (7)  form  a  elosed  set  of  equations  which  can  be  solved  to  obtain 
eaeh  component  of  the  velocity  perturbations  in  terms  of  a  single  set  of  coefficients  a„.  The 
vorticity  modes  are  dependent  on  x,  but  the  x  dependence  is  small  for  self  similar  flows  (such  as 
2D  wakes, axisymetrie  wakes  and  mixing  layers)  whieh  satisfy  the  condition  that 

(8) 


where  fTand  A  are  constants  and  L  is  the  wake  width. 

The  turbulent  shear  stresses  ean  be  obtained  by  taking  the  expected  value  of  the  velocity 
components  and  noting  that  if  the  mode  amplitudes  are  uncorrelated  then  the  velocities  arc 
statistically  stationary  in  time  and  locally  homogeneous  in  the  x  and  z  directions.  We  also  need  to 
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integrate  over  all  frequencies  and  spanwise  wavenumbers  to  obtain  the  fluctuations  at  a  point. 
Applying  these  steps  leads  to  a  description  of  the  turbulent  stress  in  the  upwash  direction  as 

(9) 


00  00  OO 

ojy)=\  J  X  A„m0)Vl”\y)d(od. p  r^Hv^OOl2 

n= 1 


similarly  for  the  other  two  components  of  the  flow  we  have, 

P2  f  x  du  V 


o»u  +  =  J  J  £  A„  (0),  p)\  U(n) (>’)  + 

-00  \ 


U  dy  J 


V(n\y)82{(0,x,y)  \daxip 


.2 

d(Ky) 


(10) 


the  shear  stress  is  then 


O',...  —  — 


(  j  j  \  °°  00  ^  ,  j  2 

JJ-Z-  J  JX  -jrAi^P)v(,,\y)^(S{(i),x,y))d6)dp 

v  ;»=1  A,n 


(ID 


U  dy 

To  a  good  approximation  Re(^)  =|  8  |2,  so  we  can  simplify  (10)  and  specify  the  tke  as 

jc  dU  ^ 


(12) 


KOO-rJ  J  I  An{(0,P)(U{n){y)  +  VM{y))d(0dP -i 

^  n=\  L 


U  dy  ) 


o..v 


If  we  ean  solve  this  equation  for  the  unknown  coefficients  A„  then  we  ean  use  the  result  to 
obtain  eaeh  of  the  stresses  and  the  two  point  correlation  of  the  velocity  fluctuations.  For 
example  the  two  point  correlation  of  the  upwash  velocity  is  given  as 

(13) 


/U^/)=  J  J  X  An(0), pWn\y)v(n) (y')doxip 

_  »=i 


The  analysis  given  above  describes  how  all  the  turbulent  flow  quantities  of  interest  ean  be 
specified  in  terms  of  a  single  set  of  mode  amplitudes  A„  .  To  bring  the  system  to  elosure  it  is 
neeessary  to  invert  any  one  of  the  above  equations  to  obtain  the  mode  amplitudes.  However  the 
mode  amplitudes  are  also  dependent  on  the  frequency  and  the  spanwise  wavenumber  and  so  it  is 
also  neeessary  to  assume  so  form  for  the  energy  distribution  of  the  turbulenee. 

1.3  Approximate  Solution  for  the  Continuous  Modes 


In  many  flows,  especially  the  wake  flow  which  is  of  interest  here,  wc  can  make  the 
approximations  that  the  Reynolds  number  is  very  large  and  that  ?2JU  »  d2U  fdy2  and  so  the  last 
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two  terms  in  equation  (1)  are  small  in  comparison  with  the  other  terms  giving  the  approximate 
equation 

(14) 

J2vin) 

(o-anU){^-Xyn))^  0 
dy 


The  solution  to  this  equation  is  obtained  directly  using 

i(o-anu(y)Wn)(y)=o 


(15) 


^2  («) 

(~-p — Xnvin))«q(n) 


from  which  it  follows  that 


(16) 


^"')(y)  =  (t5(y->’n)  vin\y)  =  I  Giy\y"rf'Xy')dy'=\iG(y  ly„) 


where  y„  is  the  solution  to  the  equation  Qf=anU(yn)  ,  G(ylyn)  is  the  Greens  function  for  Poissons 
equation  which  satisfies  the  boundary  conditions  of  the  flow,  and  (j.  is  a  sealing  factor  determined 
by  (2)  and  is 

_  (17) 

I  7 


J  IGCylyJPdy 


It  is  then  relatively  simple  to  complete  the  integrals  needed  in  equations  (9)-(13)  using  these 
modes  and  the  approach  given  by  Glegg  et  al  2009. 


1.4  Determination  of  Mode  Amplitudes 

One  of  the  important  features  of  the  turbulent  velocity  modes  in  the  wake  is  that  they  can  also  be 
used  to  define  the  turbulent  kinetic  energy  the  wake,  as  given  by  equation  (12).  The  turbulent 
kinetic  energy  is  an  output  from  RANS  calculations  and  so,  in  principle,  we  can  estimate  A„  if 
lfn>  and  V*”1  are  known.  The  approach  used  here  is  the  solve  the  Orr  Sommerfeld  equations  for 
the  wake  flow  to  obtain  the  modes,  and  use  these  as  a  basis  for  solving  of  the  mode  coefficients. 
However,  some  additional  modeling  is  required  to  give  the  spectral  content  of  the  modes.  To 
achieve  this  we  assume  that  the  coefficients  An((0,p)  can  be  written  as  the  product  of  a  universal 
energy  spectrum  E(XL)  and  a  mode  amplitude  in  the  form 

_  (18) 

An((0,P)  =  AnEanL)  A„=V«;+/*: 

The  universal  energy  spectrum  is  modeled  by  a  von  Karman  energy  spcctnim  (sec  Glegg  et  al 
(2008))  and  is  consistent  with  the  energy  cascade  of  turbulent  scales.  It  depends  on  a  lengthscalc 
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L  which  also  has  to  be  modeled  or  calculated  from  the  dissipation  in  the  RANS  calculation 
(Glegg  et  al  (2008)).  If  the  turbulent  kinetic  energy  is  known  at  ym  locations  in  the  flow  and  the 
mode-shapes  are  determined  by  the  solution  to  the  Orr  Sommerfeld  equation  then  we  can  use  a 
collocation  scheme  to  determine  the  mode  amplitudes  by  solving  the  equation 

(19) 

K(y.)=±  \<L.  &.=4i  /  Ea,L)Um(y.,o>,P)f  dPdm 

n=l  *■ 

Once  these  coefficients  have  been  determined  it  is  possible  to  predict  all  the  two  point  statistics 
of  the  flow  and  the  turbulent  velocity  spectra. 

1.5  Numerical  Verification  of  Results 

This  approach  has  been  used  to  analyze  the  turbulent  flow  in  a  blade  wake  based  on  the 
measurements  of  Devenport  et  al  (2001).  These  give  the  details  of  all  three  velocity  components, 
their  two  point  statistics  and  spectra  as  a  function  of  y  and  x  for  the  wake  behind  a  NACA  0012 
airfoil  at  zero  angle  of  attack  at  a  Reynolds  number  of  3.28x10s  .  The  turbulent  kinetic  energy 
(tke)  distributions  are  given  at  multiple  locations  downstream  of  the  trailing  edge  and  these  are 
found  to  be  self  similar  and  to  collapse  onto  a  single  line  as  shown  in  figure  2.  Using  the 
approach  given  above  the  mode  amplitudes  are  obtained  from  the  distribution  of  tke  and  are 
shown  in  figure  3.  To  verify  that  the  approach  is  valid  we  have  used  these  mode  amplitudes  to 
calculate  the  distributions  of  turbulent  stresses  Ovv  and  ouu,  and  the  two  point  correlation  function 
RJy,Ay)  of  the  velocity  in  the  v  direction  at  different  reference  locations  y  across  the  wake.  The 
results  are  shown  in  figures  4-6  and  are  successfully  compared  to  experimental  results.  Figure  6 
is  a  good  test  of  the  approach  because  the  correlation  function  is  not  symmetric  for  reference 
locations  that  are  displaced  from  the  wake  centerline. 

In  Figures  7  and  8  the  estimated  spectra  of  the  v  and  u  velocity  components  arc  given,  and  show 
reasonable  agreement  with  the  measurements. 


1.6  Conclusions 

In  this  study  we  have  developed  a  method  for  extracting  two  point  statistics  from  the  distribution 
of  turbulent  kinetic  energy  in  self  similar  flows.  The  approach  is  to  develop  a  modal  description 
of  the  turbulent  velocity  fluctuations  based  on  the  solution  to  the  Orr  Sommerfeld  equation  and 
use  this  to  model  the  distribution  of  turbulent  kinetic  energy  within  the  flow  in  terms  of  a  set  of 
unkown  coefficients.  If  the  distribution  of  the  tke  is  available,  cither  from  measurements  or 
RANS  calculations,  then  the  model  can  be  inverted  to  obtain  the  unknown  coefficients. 

This  approach  has  been  applied  to  the  flow  in  a  turbulent  wake  based  on  the  measurements  of 
Devenport  et  al  2001  with  very  encouraging  results.  The  model  has  successfully  predicted  the 
turbulent  stresses  and  the  two  point  correlation  function  of  the  upwash  velocity  fluctuations. 
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Part  lb:Turbulence  Modeling  for  Rotor  Stator  Interaction  Noise 

1.7  Introduction 

Turbomachinery  broadband  noise  is  often  dominated  by  the  turbulence  in  the  rotor  wakes 
interacting  with  downstream  stator  vanes.  The  accepted  approach  for  analyzing  this  problem  is  to 
model  the  turbulence  incident  on  the  stator  vanes  by  a  locally  homogeneous  turbulent  flow 
whose  spectrum  can  be  defined  by  a  von  Karman  or  Liepmann  model,  with  a  specified 
turbulence  intensity  and  length  scale.  This  is  a  very  simplistic  model  of  the  complex  turbulent 
flow  downstream  of  a  rotor,  which  is  dominated  by  rotating  wake  flows,  secondary  flows  and 
end  wall  effects.  In  an  earlier  paper  Glegg  and  Devenport  (2001(b))  proposed  a  more  general 
approach  to  turbulence  modeling  in  complex  flows  based  on  proper  orthogonal  decomposition  of 
the  unsteady  flow  incident  on  the  stator.  The  limitation  of  this  approach  was  the  determination  of 
the  proper  orthogonal  modes.  In  this  paper  we  will  revisit  this  problem  and  show  how  the  modes 
of  the  flow  can  be  obtained  from  the  solutions  to  the  linearized  Navier  Stokes  equations,  or  in  the 
case  of  quasi  two  dimensional  flows,  the  Orr  Sommerfeld  equations.  We  will  illustrate  the 
procedure  by  showing  results  for  a  simple  idealized  model  of  an  oblique  blade  wake  interaction. 

1.8  Sound  Radiation  from  an  Oblique  Blade  Wake  Interaction 

An  oblique  blade  wake  interaction  occurs  when  the  wake  from  an  upstream  blade  is  cut  by  a 
downstream  blade  at  an  angle,  as  is  typical  of  a  rotor  wake  stator  blade  interaction.  The  wake  is 
modeled  with  a  two  dimensional  mean  flow  profile  which  does  not  vary  along  the  wake  axis,  and 
the  cutting  angle  is  defined  as  the  angle  between  the  spanwise  coordinate  of  the  blade  and  the 
axis  of  the  wake  (see  Figure  1).  The  wake  coordinates  are  specified  as  (x,y,z)  and  the  mean  flow 
velocity  in  the  wake  is  U(y)  in  the  direction  of  the  x  coordinate.  The  blade  is  defined  in  terms  of 
coordinates  ( xh,yh,zh )  where  xb  is  the  chordwise  direction  and  yh  is  the  spanwise  direction.  The 
acoustic  field  from  the  blade  wake  interaction  can  be  calculated  from  the  unsteady  loading  on  the 
blade  surface,  which,  in  ducted  fan  applications,  can  be  coupled  with  the  duct  modes  using  the 
approach  given  by  Glegg  and  Walker  (1999).  Here  we  will  simplify  the  problem  by  considering 
the  radiation  from  the  blade  in  the  absence  of  the  duct.  This  gives  an  idealization  of  a  rotor  wake 
stator  interaction  which  is  suitable  for  identifying  the  principle  features  of  turbulent  wake  flow 
that  cause  broadband  fan  noise  from  rotor  stator  interactions,  without  the  complexity  introduced 
by  the  blade  response  function  for  a  high  solidity  stator  in  a  duct. 

Amiet  (1974)  gives  the  solution  to  the  general  problem  of  an  unsteady  flow  incident  on  a  blade 
in  a  uniform  mean  flow.  The  acoustic  field  is  defined  in  terms  of  the  unsteady  loading  on  the 
blade  surface,  which  is  determined  by  the  blade  response  to  a  harmonic  gust.  If  the  incoming 
gust  has  an  upwash  component  on  zh=0  of  the  form 

-ify)t+iY„xh+iv\h 

we  r°h  h 


then  Amiet’ s  blade  response  function  gives  the  unsteady  loading  on  the  blade  as 
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F(xt,ytJ)=wHxt.Y..vy-m'** 

For  the  oblique  blade  wake  interaction  problem  we  will  model  the  unsteady  flow  in  the  wake  by 
the  solution  to  the  Orr  Sommerfeld  equation  for  the  unsteady  flow  in  the  wake.  This  gives  a 
modal  solution  for  the  unsteady  flow  in  the  wake,  which  is  of  the  form 

u(w,0=|;  J  J  aH(a),pWm}(y,a),P)e^^u^zdpdto 

where  an(x)  are  the  eigenvalues  of  the  solution  and  can  be  assumed  to  be  slowly  varying  in  x,  and 
uw  are  the  eigenvectors  or  modes  of  the  solution.  This  expansion  of  the  velocity  has  some 
important  features  that  follow  from  the  statistical  nature  of  the  turbulent  flow.  We  can  assume 
that  the  turbulent  velocity  fluctuations  are  stationary  in  time  and  homogeneous  in  the  ^  direction. 
Within  the  limit  that  an(x)  is  slowly  varying  we  can  also  assume  that  the  turbulent  velocity 
fluctuations  are  homogeneous  in  the  direction  of  the  flow.  It  follows  that  the  coefficients 
an((0,fi)  are  uncorrelated  for  different  values  of  a?, wand  /Jand  consequently,  that  the  uncorrelated 
modes  for  the  wake  flow  are 

We  can  then  proceed  to  analyze  the  oblique  blade  wake  interaction  problem  in  exactly  the  same 
way  as  was  done  by  Glegg  and  Devenport  (2001(a)),  defining  the  upwash  from  the  mode 
velocity  and  the  cutting  angle  of  the  wake,  and  calculating  the  unsteady  load  and  the  acoustic 
radiation  for  each  mode  independently  to  give  the  acoustic  pressure  at  an  observer  located  at  x  as 
p<n>(\,(d).  Since  each  mode  is  uncorrelated  the  spectral  density  of  the  acoustic  pressure  is 

VX’®)  =  L  I  An((0,P)\ Pin\x,6»\2  d(i  An((0,fi)  =  ^Ex\\an((0,[i)\2] 

n~\  —  * 

The  important  feature  of  this  result  is  that  the  acoustic  field  for  each  mode  can  be  readily 
calculated,  and  the  characteristics  of  the  acoustic  field  can  be  related  to  the  features  of  each 
mode,  but  it  remains  to  determine  the  mode  amplitudes. 

1.9  Acoustic  Predictions 

Using  the  results  of  Part  1(a)  it  is  possible  to  carry  out  acoustic  predictions  for  the  sound 
radiation  from  a  blade  wake  interaction  at  different  cutting  angles.  In  figure  9  we  show  some 
preliminary  results  of  the  effect  of  cutting  angle  on  the  far  field  noise.  We  note  that  as  the  cutting 
angle  is  reduced  the  radiated  sound  level  is  increased  significantly.  This  is  caused  by  two 
different  effects.  First,  at  smaller  cutting  angles  the  wetted  surface  area  of  the  blade  wake 
interaction  is  larger,  and  secondly  the  upwash  is  impacted  by  different  components  of  the 
turbulent  velocity  as  the  cutting  angle  is  changed.  For  a  perpendicular  blade  wake  interaction 
(90  deg  cutting  angle)  the  upwash  on  the  blade  is  determined  by  the  z  component  of  the 
turbulence  in  the  wake.  For  a  45  deg  cutting  angle  the  upwash  is  determined  equally  by  the  z  and 
the  y  component  of  the  upwash.  To  accurately  assess  the  importance  of  the  cutting  angle  on  the 
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radiated  noise  it  is  necessary  to  determine  the  relative  amplitude  and  wavenumber  spectra  of 
these  two  different  velocity  components.  The  modeling  approach  described  here  allows  this  to  be 
done  by  utilizing  the  solution  to  the  Orr  Sommerfeld  equation  for  the  unsteady  flow  in  the  wake. 
This  is  a  far  more  robust  approach  than  one  that  assumes  a  homogeneous  turbulent  flow  model. 

1.10  Conclusions 

This  paper  describes  a  new  approach  for  describing  inflow  turbulence  for  leading  edge  noise 
problems.  The  unsteady  inflow  is  modeled  by  the  Navier  Stokes  equations  linearized  relative  to 
the  mean  flow  velocity.  The  resulting  solutions  give  the  relative  magnitudes  of  all  velocity 
components  throughout  the  flow,  and  this  provides  a  model  of  the  local  tke  in  the  flow  that  can 
be  matched  to  the  tke  calculated  by  RANS  models.  The  resulting  unsteady  flow  model  can  then 
be  used  to  carry  out  acoustic  predictions  of  leading  edge  noise.  This  approach  has  been  shown  to 
successful  match  measurements  of  two  point  correlations  and  spectra  of  the  turbulent  flow  in  the 
wake  of  a  NACA  0012  airfoil,  and  the  application  to  leading  edge  noise  has  been  demonstrated. 
The  extension  of  this  method  to  ducted  flows  and  the  wake  of  a  rotor  is  also  possible  based  on 
the  model  of  the  linearized  flow  in  a  circular  duct  given  by  Atassi  (2005). 
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Figure  1:  the  coordinate  system  for  a  simple  wake  incident  on  an  isolated  blade 
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Turbulent  Kinetic  Energy 


Figure  2:  The  collapse  of  the  tke  in  a  turbulent  wake  behind  a  NACA  0012  airfoil  as  a  function  of 
displacement  from  the  centerline.  The  data  is  self  similar  an  collapses  on  the  wake  half  w  idth  L.  The  mean 
line  is  the  data  fit  obtained  from  the  modal  solution. 
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Figure  3:  The  mode  amplitudes  calculated  from  the  distribution  of  kinetic  energy 
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Figure  4:  The  estimated  distribution  of  u:  across  the  wake. 
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Figure  5:  The  estimated  distribution  of  v2  across  the  wake. 
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Figure  6:  The  two  point  correlation  function  Rn(y,Ay)  for  different  reference  locations  y  across  the  wake 
(a)  y=0.5L.  (b)  y=L  (e)  y=2L 
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Figure  7:  The  spectra  G„  at  two  different  locations  across  the  wake,  (a)  y=L  (b)  y=2L 
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Figure  8:  The  spectra  G„„  at  two  different  locations  across  the  wake,  (a)  y=L  (b)  y=2L 
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Figure  9:  Predicted  aeoustie  spectra  for  leading  edge  noise  from  an  oblique  blade  wake  interaction.  Spp  is 
given  in  dB  with  an  arbitrary  reference.  Flow  Mach  number  is  0.5. 
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Part  2:  The  Tailored  Greens  Function  for  a  Step 


2.1  Introduction 

In  a  recent  article  Ji  and  Wang(2010)  described  detailed  Large  Eddy  Simulations  (LES)  of  boundary  layer 
flows  over  forward  and  backward  facing  steps.  One  of  the  objectives  of  this  study  was  the  calculation  of 
the  sound  radiated  by  the  interaction  of  the  flow  with  the  step.  To  carry  out  their  calculations  Ji  and  Wang 
made  use  of  Howe’s(1989)  low  frequency  approximation  for  the  Greens  function  tailored  to  the  step 
geometry.  This  Greens  function  gives  a  far  field  directionality  that  is  equivalent  to  a  streamwise  dipole 
and  predicts  there  will  be  a  null  in  the  acoustic  far  field  for  an  observer  directly  above  the  step.  In  contrast 
experimental  measurements  (Farabce  and  Zoccola(1998),  Catlett(2010))  indicate  that  this  is  not  the  case 
and  the  directionality  of  the  sound  is  relatively  weak  at  low  frequencies.  Catlett(2010)  also  observed  that 
there  was  an  interference  effect  in  the  spectrum  at  high  frequencies  that  could  not  be  predicted  by  the  low 
frequency  approximation.  In  this  paper  the  far  field  approximation  for  the  Greens  function  tailored  to  the 
step  geometry  will  be  derived  by  applying  the  Weiner  Hopf  method.  The  result  shows  that  the  sound  in 
the  acoustic  far  field  is  only  equivalent  to  a  dipole  at  very  low  frequencies,  and  it  has  very  different 
characteristics  at  higher  frequencies.  The  far  field  sound  spectrum  can  also  be  estimated  and  it  is  shown 
that  the  interference  effect  observed  by  Catlett(2010)  is  correctly  indentified. 

2.2  The  Greens  Function 


The  geometry  for  a  step  of  height  b  is  shown  in  Figure  1  in  terms  of  the  Cartesian  coordinates 
The  origin  of  the  coordinate  system  is  at  the  bottom  of  the  step  and  the  surface  is  defined  by  y:=0,  y/>0 
and  V2=b,y/<0.  To  calculate  the  far  field  Greens  function  we  will  use  the  adjoint  method,  which  utilizes 
the  principle  of  reciprocity.  In  this  approach  the  source  is  placed  at  a  location  in  the  acoustic  far  field  and 
the  observer  is  placed  in  the  vicinity  of  the  step.  The  advantage  of  using  this  method  is  that  the  adjoint 
Greens  function  can  be  obtained  by  calculating  the  field  scattered  by  the  step  due  to  sound  from  the  far 
field  source,  and,  local  to  the  step,  the  incident  field  is  well  approximated  by  a  plane  wave.  The  Greens 
function  for  the  source  close  to  the  step  and  the  far  field  observer  is  then  simply  the  adjoint  Greens 
function  with  the  source  and  observer  positions  reversed. 


In  the  adjoint  problem  the  source  is  placed  at  jc,-  and  the  incident  wave  is  given  by 

ex p ( ikr- ikxyJM nr  where  r=|x|  is  the  distance  of  the  source  from  the  bottom  of  the  step  and  k  is  the 
acoustic  wavenumber  at  the  frequency  CO.  The  adjoint  Greens  function  which  satisfies  the  non  penetration 
boundary  conditions  on  y2=0  is  then  given  by 

(1) 


^  ikr-ikxy\y  Ir-ikx  |  V|  /r 

GJy,  I  x,)  = - - - cos (kx2y2  /  r)  r  »l  y  I 

2nr 


To  satisfy  the  boundary  conditions  on  the  step  we  need  to  add  solutions  that  correct  for  the  non 
penetration  boundary  condition  on  the  surfaces  y2=b,yi<0  and  0 <y2<bf  yi-0 .  We  will  treat  these  as  the 
sum  of  two  separate  corrections  Gj  and  G>  First  wc  define  G/  as  the  solution  that  corrects  G0  to  match  the 
boundary  conditions  on  the  upper  surface  of  the  step,  so  it  has  the  boundary  conditions  (see  figure  2) 

(2) 
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Secondly  we  define  G’  as  the  solution  that  corrects  G„+Gi  to  match  the  boundary  conditions  on  the  face 
of  the  step,  so  it  has  the  boundary  conditions  (sec  figure  3) 
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The  solution  for  Gj  is  obtained  from  the  result  given  by  Noble(1958)  for  the  scattering  of  sound  by  a  semi 
infinite  two  dimensional  duct.  If  we  represents  the  Greens  function  G0  as  the  sum  of  two  waves  incident 
from  the  angles  0and  - 0 ,  where  cos6cos<p=x j/r,  sin{fcos^=x/r  and  sin0=v yh  Then  we  can  write 

(4) 

ikvr—ikxy\ylr 
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and  use  the  result  given  by  Noble  (1958)  for  scattering  of  plane  waves  by  two  semi  infinite  planes.  The 
result  is  given  by  Noble’s  equation  3.25  as 


G,  =  —  f  SAa)E(yy2)e~^da  E(yy2)  = 

IT  * 


f-cosh(yv,  )e 
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(7) 
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The  choice  of  branch  cut  for  the  square  roots  is  important  and  we  choose  the  branches  where 
Re((Cf^rrt)1  ^>0.  The  subscripts  ±  represent  functions  which  are  analytic  for  Im(a)  greater  than  or  less 
than  zero  respectively,  and  specifics  of  L+  and  L  will  be  discussed  later. 

The  next  step  is  to  determine  the  second  correction  which  ensures  the  velocity  on  the  step  face  is  zero.  We 
note  that  the  sum  of  G0  and  G}  result  in  waves  which  propagate  from  right  to  left  along  the  duct  defined 
by  0 <y2<by  yi<()  (  see  Figure  3).  To  cancel  the  velocity  at^/=0  wc  introduce  a  wave  which  propagates 
from  left  to  right  along  the  duct  so  that  the  boundary  condition  is  satisfied.  The  new  wave  will  also  be 
transmitted  across  the  duct  termination  and  propagate  out  of  the  duct  and  so  will  affect  the  region  of 
interest.  This  problem  is  also  discussed  by  Noble  and  for  a  wave  Bno,os(rm(y2-b)/b)cxv(iKn)'i)i  (where 
Kn=( k02- (nn/b)2) !  \  and  has  a  positive  imaginary  part)  propagating  along  the  duct.  The  acoustic  field  is 
given  by  Noble  equation  (3.34)  as 

(8) 
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To  obtain  the  unknown  coefficients  B„  we  evaluate  the  field  in  the  ducted  region.  Noble  shows  that  by 
evaluating  equation  (4)  we  have 

(10) 
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=  X  cm  cos (mn(y2  -  b)  /  b)e"K"y' 
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where  Cm  represents  the  residue  at  the  poles  of  the  integrand  of  (5)  which  lie  in  the  upper  half  of  the 
complex  plane  and  are  given  in  the  appendix.  Similarly  we  can  evaluate  (8)  in  the  dueted  region.  The 
solution  includes  the  waves  propagating  from  left  to  right  and  the  waves  which  are  reflected  from  the  duct 
termination  at  yj=0,  and  ean  be  expressed  as 

(ID 

[C,L,<.,<.  =  X  I  B.(S„e*-»+Rm,e-“~-)cosimUy,-b)lb) 

m= 0  n=0 

where  the  reflection  coefficients  Rmn  are  also  given  in  the  appendix.  Using  (10)  and  (1 1)  in  the  last  of  the 
boundary  conditions  (3)  and  equating  terms  gives 

(12) 

I  BjSmn-Rnm)  =  Cm 

n~  0 

If  this  series  is  truncated  to  a  finite  number  of  terms  then  this  equation  can  be  solved  for  the  coefficients 
Bn ,  giving  the  complete  solution  for  G\ 

2.3  Approximate  Solutions 


At  very  low  frequencies  we  ean  show'  that  L+(k0)^l  and  Bn~~Co/2=A  and.  Also  the  series  expansion  for  G: 
is  dominated  by  the  n=0  term  and  so  for  v:>b  we  ean  approximate 

(13) 
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Applying  the  same  approximations  to  G/  in  the  region  }';>h  ,  and  utilizing  equation  (7)  gives 


(14) 
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At  moderate  frequencies  where  k„b=l  the  situation  is  quite  different.  In  this  ease  we  find  from  numerical 
evaluation  that  B„~C„  and  that  Bn«B„ ,  but  we  can  not  approximate  the  other  terms.  Then,  for  )’2>b  we 
obtain 
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and  wc  note  that  the  scaling  on  k„  and  directionality  associated  with  these  functions  is  quite  different  from 
the  low  frequency  approximation. 

2.5  Discussion 

The  sound  radiation  from  turbulent  flow  near  rigid  surfaces  can  be  evaluated  by  using  Lighthills  acoustic 
analogy.  The  sound  sources  arc  specified  by  Lighthills  stress  tensor  and  the  acoustic  field  is  given  by 
(Goldstein  (1976)) 

(17) 

P'c]{xr(D)=\  Ttj ( V, , Q) ) dV 

v  v\,  o)j 

where  G  is  a  tailored  Greens  function  that  satisfies  the  non  penetration  boundary  condition  on  all  the 
surfaces.  The  radiation  efficiency  of  the  sound  sources  is  therefore  given  by  d2 G/dy,dyj , which,  for  flow 
over  a  step  is  given  by  the  Greens  function  defined  above.  If  we  evaluate  fr GJdytdyj  it  is  found  that  the 
aeoustic  radiation  scales  as  k“  and  this  leads  to  a  souree  of  quadrupole  order.  In  contrast  if  we  evaluate 
oGi/dyi  we  find  that  this  simple  scaling  is  not  present.  If  we  model  the  Lighthills  stress  tensor  for  a 
turbulent  eddy  convetced  in  the  V/  direction  as 

(18) 

Tlj(yi>Q))  =  fijel°»',u5(y2-h)f 

where  h>b  then,  and  t  is  a  length  scale,  then  using  the  low  frequency  approximations  we  obtain 
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where 


(21) 
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x„  =  V(«/t/)2-^2 


At  very  low  Maeh  numbers  where  (tX/U»k„  we  ean  approximate  y„~(o/U  and  it  follows  that 

(22) 

p'cl(x,,(o)~{-2iAknbcosO)(o)f/U)L+(co/U)(-fu-ifr_+f22) 


First  we  note  that  this  result  has  the  eharaeteristies  of  a  streamwise  dipole,  sealing  with  kQ  and  having  a 
directionality  of  cosft  The  sealing  on  flow  speed  ean  be  obtained  by  noting  that  scales  as  pU2  and  k0 
scales  as  M(dU  where  M=U/c0 ,  and  ca  is  the  speed  of  sound.  Then  at  constant  Strouhal  number  (cob/U) 
the  acoustic  pressure  scales  as  if  and  the  acoustic  power  output  scales  as  if. 

For  the  higher  frequencies  the  sealing  is  quite  different.  Equations  (15)  and  (16)  show  that  the  sealing  on 
frequency  depends  on  s\n(kobs\n0)  and  L+(k0cosG),  so  the  field  no  longer  has  the  eharaeteristies  of  a 
dipole. 

The  complete  eharaeteristies  of  the  directionality  of  the  sound  generated  by  flow  over  a  step  can  be 
calculated  directly  from 

(23) 

p'c2o(x„0))  = 

f  2/A  sin  OsinQ^b  sin  /(/  +  £„)  ,y  +A:J(ft)/(7  +  A:n)^.(-vn)] 

\(\  +  cosd)L+(k0cos6)(o)  /  U  -  kacos6)  (o)IU  +  Kn)  [ 

x  (f/yo)L+(ct)  /  /  U)2fu  -  i((0  /U)Y„ fl2  +  ylt2)e~r^ 


This  result  has  been  evaluated  to  show  the  directionality  of  the  sound  field  at  different  non  dimensional 
frequencies  kob-0.01 ,0.1 ,1 ,  For  a  backward  facing  step  in  a  low  Mach  number  flow  IM-0.02)  are  shown 
in  figure  4.  At  low  frequencies  the  directionality  has  a  minimum  at  0=90°  and  is  similar  to  a  dipole 
orientated  in  the  direction  of  the  flow.  At  high  frequencies  the  directionality  is  impacted  by  interference 
effeets  and  is  almost  omnidirectional.  Similar  results  are  shown  for  a  forward  faeing  step  in  Figure  5, 
(obtained  from  equation  (23)  w  ith  U  replaced  by  -U)  and  it  is  hard  to  distinguish  the  difference.  However 
at  higher  Maeh  numbers  ( M-0.2 )  the  directionality  is  different  as  shown  in  figures  6  and  7,  and  the  effeet 
of  flow  direction  is  more  marked  w  ith  nulls  pointing  tow  ards  the  direction  of  the  flow. 

Figure  8  show  s  the  frequency  dependence  of  the  field  at  three  different  observer  angles  for  a  forward 
facing  step  at  a  Mach  number  of  0.2.  It  is  impotant  to  note  the  interference  dip  that  occurs  when 
knbQQ$9=n  and  the  lack  of  collapse  of  the  aeoustic  efficiency  onto  a  single  eurve.  This  indicates  that 
sealing  of  experimental  data  on  a  single  frequency  parameter  could  be  unsuccessful. 

In  Figure  9  show  the  spectrum  of  the  far  field  sound  is  shown  assuming  a  source  spectrum  w  hich  seales  as 
(kjh)'x  x  for  an  1  1  mm  step  and  observer  positions  that  correspond  to  those  measured  by  Catlett(2010).  An 
arbitrary  scaling  factor  has  been  added  to  align  the  plots  so  this  is  not  an  absolute  prediction.  However  he 
spectra  show  the  interference  dip  and  resonances  at  frequencies  above  10kHz  which  were  observed 
experimentally  by  Catlett  (2010)  and  there  is  some  consistency  between  the  measured  levels  at  different 
angles  at  low  frequencies. 
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2.6  Conclusion 


The  results  given  above  discuss  the  sound  radiation  from  turbulent  flow  over  a  step.  It  is  concluded  that 
the  sound  radiation  is  a  consequence  of  scattering  mechanism  that  can  be  modeled  from  solutions 
presented  by  Noble  for  scattering  by  parallel  semi  infinite  plates.  At  very  low  frequencies  the  far  field 
sound  has  the  characteristics  of  a  streamwise  dipole  and  the  sound  from  turbulent  flow  in  the  vicinity  of 
the  step  scales  with  the  sixth  power  of  the  flow  velocity.  At  higher  frequencies  the  there  interference 
effects  become  important  and  the  spectrum  exhibits  an  dip  at  angles  where  kohsin&^nn:.  The  directionality 
and  spectral  characteristics  are  similar  to  those  measured  experimentally  by  Catlett(2010). 
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Appendix 


The  coefficients  Cm  are  obtained  from  (5)  as 

_  2As\n(kobs\n0) 

fc>stn0L+(fc„cos0)L_(fc„) 

c  _  _  -2A  sin  flsin(fc„/?sin  0)Lt  (*•„  )(mnf _  ^  >  Q 

Kmb(\+cos0)(Km-kocose)(Km-ko)b2L+(kocosG) 

Similarly  we  obtain  for  the  reflection  coefficient: 

R  .  t L )U~- + *»  )L-(-k~«  ) 

where  sm-l  for  m>0  and  eo~0.5. 


The  Split  Function 

Following  Noble  pi 04  we  can  define  the  split  function 
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and  L.(a)=L+(-a). 
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Figure  1 :  The  eoordinate  system  used  to  define  the  step  and  the  adjoint  problem 
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far  field  source 


Image  source 


Figure  2:  Schematic  of  the  first  decomposition  which  satisfies  the  boundary  conditions  on  V2=k,yi<0  and 
y2~0,yj>0. 
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Figure  3:  Schematic  of  the  second  problem  that  satisfies  the  boundary  conditions  on  the  face  of  the  step. 
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Figure  4:  The  Directionality  of  sound  from  sources  near  a  backward  facing  step  at  a  Mach  number  of 
0.02.  Levels  are  normalized  by  kQb  and  50dB  has  been  added  to  ensure  positive  values  on  all  plots,  -x-x- 
k„b~0.01,  -o-o-kob~01 ,  knb-l 
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Figure  5:  The  Directionality  of  sound  from  sources  near  a  forward  facing  step  at  a  Mach  number  of  0.02. 
Levels  are  normalized  by  kftb  and  50dB  has  been  added  to  ensure  positive  values  on  all  plots,  x-x- 
kob=0.01,  -o-o-kob=01 ,  kQb=l 
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Figure  6:  The  Directionality  of  sound  from  sources  near  a  backward  facing  step  at  a  Mach  number  of  0.2. 
Levels  are  normalized  by  k„b  and  50dB  has  been  added  to  ensure  positive  values  on  all  plots,  -x-x- 
kob=0.0],  -o-o-knb=0 1 ,  k„b=l 
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Figure  7:  The  Direetionality  of  sound  from  sources  near  a  forward  facing  step  at  a  Mach  number  of  0.2. 
Levels  are  normalized  by  kab  and  50dB  has  been  added  to  ensure  positive  values  on  all  plots,  -x-x- 
kob=0.01,  - o-o-kj)=01 ,  k0b=l 
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Spectrum  Level 


Figure  8:  The  Greens  function  as  a  function  of  kah  for  a  forward  facing  step  at  a  Mach  number  of  0.2  at 
three  different  angles  6^56. 5°, 82.5°, 1 06°,  and  128.5° 
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Far  Field  from  a  Forward  Step  of  h  =  11.7  mm  and  U  =  60  m/s 


Figure  9:  Comparison  of  measurements  by  Catlett(2010)  and  the  predictions  based  on  the  Greens  function 
developed  in  this  study. 
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Part  3:  Predictions  of  Surface  Pressure  Spectra  for  a  Turbulent  Boundary 
Layer 

3.1  Introduction 

This  part  of  the  report  will  discuss  a  method  for  extracting  velocity  and  surface  pressure  spectra 
associated  with  a  high  Reynolds  number  turbulent  boundary  layer.  The  approach  is  to  develop  a 
description  of  the  turbulent  velocity  fluctuations  based  on  the  solutions  to  the  Orr  Sommerfcld 
and  Squire  equations  and  to  use  these  to  model  the  distribution  of  turbulent  stresses  within  the 
flow  in  terms  of  a  set  of  unknown  coefficients.  If  the  distribution  of  the  turbulent  stresses  is 
available,  either  from  measurements  or  RANS  calculations,  then  the  model  can  be  inverted  to 
obtain  the  unknown  coefficients  and  closure  is  obtained. 

This  approach  has  been  applied  to  the  flow  in  a  turbulent  boundary  layer  downstream  of  a 
forward  facing  step  based  on  the  measurements  of  Awasthi  et  al  2011.  The  modeling  has 
identified  terms  associated  with  large  scale  structures  in  the  flow  and  small  scale  turbulence,  and 
an  energy  budget  for  each  term  has  been  obtained  as  shown  in  Figure  7.  Another  output  of  the 
model  is  the  turbulent  velocity  spectra  and  these  are  compared  to  measured  spectra  at  three 
heights  above  the  wall  in  Figure  8.  The  interesting  result  is  that  by  using  the  viscous  equations 
with  non  linear  interactions  allowed  for,  the  spectra  arc  estimated  over  a  wide  range  of 
frequencies  without  recourse  to  a  modeled  turbulence  energy  spectrum. 


3.2  Theoretical  Background 

In  the  following  we  will  consider  methods  for  evaluating  the  unsteady  part  of  a  fully  developed 
turbulent  tlow.  We  will  limit  consideration  to  incompressible  constant  temperature  flows  and 
assume  that  the  mean  flow  is  known  from  a  solution  to  the  Reynolds  averaged  Navicr  Stokes 
equations.  The  flow  is  defined  by  its  steady  velocity  and  pressure  U  and  P  and  its  unsteady 
velocity  and  pressure  u  and  p.  The  equation  of  motion  describing  the  unsteady  flow  is  obtained 
by  subtracting  the  Reynolds  average  equation  for  the  steady  flow  from  the  full  Navier  Stokes 
equations,  to  give 

(1) 

—  +  (U.V)u  +  (u.V)U  +  —  V/j — —  Vu  =  -b  V.u  =  0 

dt  p  R 

b  =  (u.V)u  -  ((u.V)u) 

The  left  side  of  this  equation  is  linear  in  the  unsteady  flow  variables  and  the  vector  b  specifies 
the  contributions  from  the  non  linear  terms.  In  the  following  we  will  also  need  a  Poisson 
equation  for  the  pressure,  which  is  obtained  in  the  usual  way  by  taking  the  divergence  of 
equation  (1),  and  the  vortieity  equation,  which  is  obtained  fro  the  curl  of  equation  (1) 


37 


Wc  will  consider  velocity  perturbations  u  =(u,v,w)  in  a  parallel  shear  flow  defined,  in  Cartesian 
coordinates  ( x,y,z )  by  the  mean  velocity  U(y)  in  the  x  direction.  The  linearized  form  of  the 
Navier  Stokes  equations  and  the  continuity  equation  are  used  to  obtain  the  Orr  Sommerfeld 
equation  (see  for  example  Hallbaek  1996)  by  taking  the  Laplaeian  of  (1)  for  the  velocity 
component  in  the  y  direction,  subtracting  Poissons  equation  for  the  pressure,  and  ignoring  the 
non  linear  term.  However  since  the  non  linear  terms  are  important  in  a  fully  developed  turbulent 
flow  we  will  retain  them  here  to  give  an  inhomogeneous  form  of  the  Orr  Sommerfeld  equation  as 

(2) 


— (V2v)  -U"—-  —  V4v  =  -V26  +  (V.b)' 
Dt  dx  R  ' 


where  the  primes  represent  differentiation  with  respect  to  they  direction,  D,/Dt=d/dt+(Jd/dx,  and 
b=(bXlby,b2).  The  u  and  w  components  of  the  velocity  ean  be  obtained  by  considering  the  vortieity 
in  they  direction,  given  by  r]=du/dz-dw/dx,  and  the  continuity  equation  (Hallbaek  (1996)).  The 
equations  for  the  u  and  w  velocity  components  are  then  given  in  terms  of  the  normal  velocity  and 
the  normal  vortieity  as 


d~u  d~u  c?"v  dt]  d'w  d~w  d2v  dt] 

dx 2  dz1  dxdy  dz  dx2  dz2  dzdy  dx 


(3) 


To  close  this  set  of  equations  we  use  the  vortieity  equation  to  obtain  an  inhomogeneous  form  of 
Squires  equation  for  the  normal  component  of  the  vortieity  as 


M 

Dt 


V2  77  +  (/'^  =  -Vxh 

Rt  dz 


(4) 


These  homogeneous  form  of  these  equations  have  been  considered  by  many  different  authors 
(see  for  example  Drazin  and  Reid(  1 98 1 ),  Hallbaek  (1996)  for  reviews)  to  investigate  flow 
instabilities  and  boundary  layer  transition.  Solutions  have  been  sought  that  identify  the  growth  or 
deeay  of  disturbances  and  it  ean  be  shown  that,  in  the  inviseid  ease,  the  amplitude  of  the  normal 
component  of  the  velocity  grows  linearly  with  time  from  an  initial  disturbance.  Also  it  ean  be 
shown  that  a  boundary  layer  flow  will  be  unstable  if  the  mean  velocity  profile  includes  a  point  of 
inflexion.  However  in  this  study  we  are  interested  in  using  linear  theory  to  describe  a  high 
Reynolds  number  fully  developed  turbulent  flow,  which  is  characterized  as  stationary  in  time  and 
homogeneous  in  the  x  and  z  directions.  On  this  basis  we  will  seek  solutions  in  the 
frequeney/wavenumber  domain  by  considering  the  Fourier  transform  of  all  variables,  defined  as 

(5) 

.  T  R  R 

f(y,a),a,P)  =  ——jj  J  J  / \x, y,  z, t)e~,(0,~iax~'Pz dxdzdt 

\LK)  _T  _R  _R 

Taking  the  Fourier  transform  of  the  inhomogeneous  Orr  Sommerfeld  and  Squire  equations  then 

gives 

(6) 
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(aU  +  a>)(v"-  A2v)  -  a\J"v  -  2  A2  v"  +  A4v)  = 

( aU  +  Q))ri  +  pu'v-—(fj"-X1fj)  =  C 

R< 

A2  =a2  +  P2  B  =  iX2b +  aft'  +  pb[  C  =  -(Jbx  +  ab 
and  the  u  and  w  components  of  the  velocity  are  defined  as 

(7) 

.  -i(av'-pfi)  A  -i((iv'+ art) 

U  = - r -  W  = - r - 

A2  A2 


This  provides  a  closed  set  of  equations  which  may  be  solved  to  obtain  all  the  velocity 

components  of  the  flow  in  terms  of  the  unknown  functions  B  and  C .  We  w  ill  start  by  solving  for 
the  normal  velocity,  which  can  be  considered  independently,  and  then  consider  the  normal 
vorticity,  which  depends  on  the  distribution  of  normal  velocity.  With  these  two  results  wc  can 
then  obtain  the  velocity  components  parallel  to  the  mean  flow  given  by  equation  (7). 


3.3  The  Normal  Velocity 


In  order  to  solve  equation  (6)  we  need  to  specify  the  details  of  the  mean  velocity  profile  in  the 
boundary  layer.  This  has  a  discontinuity  in  its  slope  at  the  interface  between  the  viscous  sublayer 
and  the  log  layer,  and  so  we  will  treat  the  viscous  sublayer  and  the  log  layer  as  separate  regions. 
In  the  viscous  sublayer  the  non  linear  terms  on  the  right  of  equation  (6)  can  be  considered  small 
compared  to  the  viscous  terms,  and  the  velocity  profile  is  linear,  so  the  equation  (6)  reduces  to  its 
homogeneous  form  and  will  depend  on  the  boundary  conditions  on  the  surface  and  at  the 
interface  with  the  log  layer.  The  OS  equation  is  of  fourth  order  and  so  four  boundary  conditions 
will  be  required  to  specify  the  flow  in  the  viscous  sublayer.  These  are  that  v  and  v’  should  be  zero 
on  the  surface,  and  the  v  and  v’  should  be  continuous  across  the  interface  with  the  log  layer 
where  y=ys. 

In  the  log  layer  we  will  seek  a  solution  using  the  Greens  function  for  the  OS  equation  which  is 
the  solution  to 


(aU  +  a>)(G"  -  A2G)-aU"G-—(G""  -  2  A  2G"  +  A4G)  =  5(y-y, ) 

R 


(8) 


and  satisfies  the  same  boundary  conditions  as  0  at  y=ys  and  at  infinity.  Giving  the  normal 
velocity  as 

(9) 


v(y)  =  J  fl(y,  )G(y  I  y,  )d\\ 
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3.4  The  Normal  Vorticity 


The  solution  for  the  normal  vorticity  is  obtained  using  the  same  approach  as  above.  From  the 
second  equation  in  (6)  we  obtain 


(10) 


f]  =  \  (C - fiU '(y,  )v(v, ))Gn(y I  y, )dy{ 


where  G„  is  the  Greens  function  that  satisfies  the  equation 


(a 0  +  aU)Gn(y  I  y, )  - — (G"(y  I  y, )  -  A:G„(y  I  y, ))  =  S(y  -  y, ) 

A 


(11) 


In  this  case  there  are  two  contributing  source  terms.  The  first  is  from  the  non  linear  term  C  that 
acts  as  a  source  of  normal  vorticity  and  the  second  is  from  the  interaction  term  associated  with 
the  normal  velocity.  We  can  combine  equations  (9)  and  (1 1)  for  the  interaction  term  to  give 

(12) 


nt  =-J  fl(y2)Cv(yl y2)dy2 


G,(y  l.v2)  = 


J  fiU  '(v, )Gl)(y  I  y,  )G(y,  I  y2  )dy\ 


3.5  The  Greens  Functions 

The  Greens  functions  defined  above  are  shown  in  figures  1-3  for  different  frequencies  and 
wavenumbers.  In  all  cases  the  spatial  resolution  was  set  as  dy-1/500  and  the  number  of  terms  in 
the  Fourier  series  expansion  limited  to  50.  The  boundary  layer  mean  velocity  profile  was  given 
by  U=y' 7 ,  and  the  Reynolds  number  was  5x1 0\ 

In  figure  1  the  Greens  function  for  the  OS  equation  and  the  interaction  term  are  compared  for  a 
non  dimensional  frequency  of  ry=  1  and  wavenumbers  ce=-1.285  and  /J=l.lt  is  seen  that  the 
impact  of  a  source  term  (in  this  case  the  non  linear  term  B )  at  any  location  y,  has  a  non  local 
impact  on  the  flow,  which  indicates  the  possibility  of  a  large  scale  structure.  We  note  that  at  this 
frequency  the  interaction  term  that  drives  the  normal  vorticity  is  of  a  much  lower  level  than  than 
the  Greens  function  driving  the  normal  velocity.  However  at  lower  frequencies  as  shown  in 
figure  2,  the  interaction  term  dominates. 

In  figure  3  we  show  the  Greens  function  for  Squires  equation  that  shows  how  the  normal 
vorticity  is  driven  by  the  non  linear  source  term  C.  In  this  case  we  see  that  there  is  only  a  local 
reaction,  and  the  impact  of  the  source  term  is  very  local.  This  implies  that  this  term  is  responsible 
for  the  small  scale  turbulence  in  the  flow  and  is  not  associated  with  large  scale  structures. 
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3.5  Large  scale  Structures  and  Small  Scale  Turbulence 

In  the  previous  section  it  was  shown  that  the  unsteady  part  of  the  Reynolds  average  equations 
can  be  solved  by  assuming  that  the  non  linear  terms  are  small  and  act  as  sources  of  disturbance  to 
the  larger  linear  terms.  By  evaluation  the  solutions  to  these  equations  it  was  found  that  there  are 
two  components  driving  the  flow,  one  that  has  a  broad  non  local  impact  and  another  that  only 
impacts  the  flow  locally.  In  the  following  we  will  consider  these  two  different  terms  as 
representing  the  large  scale  structures  in  the  flow  and  the  small  scale  turbulence.  The  small  scale 
turbulence  only  couples  with  the  vorticity  normal  to  the  wall  and  only  has  velocity  components 
parallel  to  the  wall.  In  contrast  the  large  scale  structures  have  velocity  and  vorticity  components 
in  all  three  directions. 


We  will  assume  that  the  flow  may  be  broken  down  into  large  scale  structures  and  small  scale 
turbulence,  and  that  the  large  scale  structures  propagate  with  their  own  phase  speed  ck  and  that 
the  non  linear  terms  associated  with  each  structure  can  be  superposed  such  that 

(13) 

K 

B{t,x,y,z)  = X  AkBk{t-xl  ck,y,z) 

*=i 


where  0 <ck<] ,Ak  represents  the  amplitude  of  the  structure  and  Bk(t,y,z)  is  the  distribution  of  the 
non  linear  terms  for  the  large  scale  structure  of  order  k.  Using  this  model  we  can  define  the 
source  terms  in  equation  (6)  as 

(14) 

K 

B(y,(o,a,p)  =  £  AkBk{y,(o,p)8{a  +  (o/ck ) 

*= i 


and  it  follows  that  the  normal  velocity  is  given  by 

v(x,r)=J  J  £  Akvk (y, <y, /i ^dcodp 

where 

R 

vk (y, co,P)  =  ]  Bk(yi,(o,p)[G(y\ y, ]a=w(B/f<  dyt 

y. 


(15) 


(16) 


Similar  expression  can  be  obtained  for  r/v  and  the  velocity  components  u  and  vv  by  utilizing 
equation  (12)  and  (7). 

Summarizing  the  results  gives  the  velocity  field  as 

(17) 

K 

u (\,t)  =  £  Akuk(t-x/ ck,y,z)  +  u s(t,x,y,z) 

k=  I 
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where  u*  is  the  veloeity  perturbation  associated  with  each  large  scale  structure  and  u,  is  the 
velocity  of  the  small  scale  turbulence,  and  we  can  define 


(18) 


uk(\,t)  =  j  j  u k(y,co,p)e-UoU-x,c'}-,p:d(0(lp 
ui(x,r)=j  j  j  us(y,o),a,P)e~im~Ulx~ip:dQ)dadP 


where 


with 


(19) 


-i(avk  -pnk) 


(20) 


A 

r\k(y,a),P)  =  -  J  B*(y,,ft),/3)[Gv(y I  y, dyx 


and  for  the  small  scale  turbulence 


r«mi 

VV  = 

-iaris 

. 

a=-a)/ck 

s 

1 

l 

f?,  =  |  C(y,  )G,((y  I  v,  )dy\ 


(21) 


An  example  of  the  velocities  associated  with  a  large  scale  structure  is  given  in  Figure  4  where  it 
is  seen  that  at  this  frequency  and  wavenumber  all  three  components  are  of  similar  magnitude.  At 
lower  frequencies  the  normal  velocity  is  relatively  smaller,  and  at  higher  frequencies  the  normal 
veloeity  dominates  the  result. 

3.6  Modeling  the  Non  Linear  Terms  and  Statistical  Averages 
3.6.1  Non  Linear  Terms 

The  non  linear  terms  in  the  above  analysis  are  unknown  and  so  we  must  resort  to  some  type  of 
modeling  to  approximately  account  for  their  characteristics.  The  non  linear  terms  may  be  both 
sources  and  sinks  of  energy  and  we  must  delineate  between  these  effects.  The  reduction  of 
energy  in  the  flow  ean  be  modeled  in  the  same  way  that  we  model  non  linear  turbulent  flow 
terms  in  the  steady  RANS  equations,  by  introducing  an  eddy  viseosity  of  the  viseous  term  on  the 
left  of  the  OS  equation  and  Squires  equation.  The  sources,  or  production  terms,  remain  on  the 
right  and  it  is  these  that  we  will  model  here. 

We  will  assume  that  the  non  linear  terms  are  dominated  by  the  small  scale  turbulence,  which 
feeds  the  large  scale  structure,  and  model  the  non  linear  term  as 
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Bk=(on'^\Gn(y,yk)/Gn(yk,yk)\2 


(22) 


where  n,  and  n2  are  modeling  coefficients,  which  are  found  to  be  n,=n2=2  from  comparison  with 
experimental  data. 

3.6 2  The  Turbulent  Kinetic  Energy 

For  the  large  scale  structure  modes  we  can  define  the  turbulent  stresses  at  the  location  y  by 
averaging  velocity  at  a  fixed  point.  Since  the  flow  is  statistically  stationary  in  time  and 
homogeneous  in  x  and  y,  different  frequencies,  wavenumbers  and  modes  must  be  uncorrelated 
and  we  obtain 

(23) 

crvv(y)  =  £[l  v(x,t) I2 1  =  j  J  X  K  1  vk(y,0),P) I2  dtocip 

*= i 

Similar  expression  can  be  defined  for  the  other  velocity  components  and  the  turbulent  kinetic 
energy  for  the  large  scale  structures  is  given  by 

(24) 

*'„(>’)  =  \  J  J  I  4,2(l  uk(y,0),P)  I2  + 1  vk(y,d),P) I2  + 1  wk(y,0),P)  I2) ckod/l 

2  ^  ^  k= i 


The  total  kinetic  energy  is  then  given  by 
*-'(y)  =  /c„(.v)  +  My) 

where  ks  is  the  kinetic  energy  of  the  small  scale  structures. 


(25) 


Similarly  the  spectrum  of  the  velocity  fluctuations  is  given  by 
Sm(CD,y)=  J  X  A2k\vk(y,co,p)\2  dp 

_oo 


(26) 


To  evaluate  these  expressions  we  need  to  accurately  integrate  over  the  wavenumber  /j.  We  can 
show  that  all  terms  in  the  integral  are  a  function  of  /?  and  so  the  integrals  can  be  carried  out  over 
the  half  range  0 </k°°.  Furthermore  to  reduce  the  computational  effort  we  can  use  logarithmic 
steps  in  wavenumber  to  extend  the  integration  over  a  wider  range  of  wavenumbers  for  a  given 
number  of  computations. 

3.63  Surface  Pressure  Spectra 

The  surface  pressure  spectrum  can  be  calculated  from  the  solution  to  Poissons  equation 
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p"-X1p  =  ipavU'-pVb 


(27) 


It  is  argued  that  the  linear  term  will  be  of  order  U/u  larger  than  the  non  linear  term  on  the  right 
side  of  this  equation,  and  since  this  is  directly  proportional  to  the  normal  velocity,  it  follows  that 
the  surface  pressure  is  determined  exclusively  by  the  large  scale  structures  in  the  flow  because 
the  small  scale  turbulence  is  parallel  to  the  mean  flow.  This  has  important  modeling  implications 
because  it  implies  that  the  convection  speed  of  the  large  scale  structures  vary  from  near  zero  to 
the  outer  stream  velocity.  However  it  does  not  relate  the  scale  of  the  large  scale  structure  to  its 
convection  speed  because  nothing  has  been  assumed  about  the  scale  of  the  structure  in  the  .v 
direction. 


To  obtain  the  surface  pressure  we  first  need  to  solve  equation  (27)  subject  to  the  boundary 
condition  that  p'=0  on  y=0.  If  the  non  linear  tern  is  ignored  in  comparison  with  the  linear  term 
then  we  obtain 

(28) 


p(y)  =  77  j  ipo.il '(y,  )v(y,  )(e"Aly  +  e~My*h  1  )d\\ 

1A  0 


Then  using  the  expansion  in  equations  (17-19)  the  surface  pressure  becomes 


(29) 


P(0  =  £aJ  J  — V\)\)vMe 


-Ay, 


e-imdyxdpd(0 


a --a)  ck 


These  results  can  also  be  used  to  calculate  the  frequency  spectrum  of  the  surface  pressure 
fluctuations.  Since  the  modes  are  uncorrelated  this  is  defined  as 


oo  <*> 

A2  j 

*  =  1  -oo 


t/'(y,K(y,)«-AVy 


Ja=-(ock 


d? 


(30) 


3.7  Comparison  with  Experimental  Results 

The  model  given  here  was  compared  to  experimental  measurements  made  in  the  stability  wind 
tunnel  at  Virginia  Tech,  as  reported  Awasthi  et  al  (201 1).  We  will  present  results  for  a  fully 
developed  turbulent  boundary  layer  with  a  Reynolds  number  of  1 .73xl0\  The  mean  velocity 
profde  and  the  turbulent  stresses  as  a  function  of  displacement  from  the  wall  are  shown  in  figure 
5.  To  compare  the  measurements  to  the  theoretical  model  we  still  need  to  define  the  unknown 
coefficients  A-  To  find  these  we  can  use  equation  (23)  for  the  normal  turbulent  stress  and  use  a 
least  squares  algorithm  to  fit  the  data  to  the  model.  The  algorithm  must  be  constrained  so  that 
Ak2>0  and  the  number  of  coefficients  must  be  less  than  the  number  of  measurements  points  to 
obtain  a  converged  result.  Figure  6  shows  the  amplitude  of  each  structure  verses  mode  order  k , 
and  it  is  seen  that  the  energy  is  distributed  across  the  modes  and  no  one  mode  dominates.  The 
strongest  contribution  is  from  mode  order  2,  but  this  is  only  a  factor  of  2  above  the  other  modes. 
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With  these  coefficients  we  ean  now  estimate  the  turbulent  kinetic  energy  (tke)  from  the  large 
scale  structures  using  equation  (24).  The  result  is  shown  in  figure  7  and  is  compared  to  the 
measured  tke.  The  difference  is  attributed  to  the  tke  from  the  small  scale  structures  and  this  is 
also  shown  in  figure  7.  The  result  indicates  that  the  energy  is  the  small  scale  structures  dominates 
the  tke  budget  by  a  factor  of  about  two. 

As  an  independent  check  on  the  validity  of  the  results  we  ean  also  calculate  the  spectrum  of  the 
normal  velocity  using  equation  (26)  and  this  is  compared  to  the  measured  results  in  figure  8  at 
three  different  locations  in  the  boundary  layer.  The  speetral  shape  appears  to  be  well  modeled  by 
this  approach,  and  is  an  interesting  result  beeause  no  speetral  modeling  or  lengthscales  have  been 
required  to  make  this  ealeulation. 

The  prediction  of  the  unsteady  flow  have  also  been  applied  to  a  boundary  layer  downstream  of  a 
forward  faeing  step  as  measured  by  Awasthi  ct  al  (201 1).  Figure  9  shows  the  boundary  layer 
profile  and  the  distributions  of  turbulent  stresses  at  nine  step  heights  downstream  of  the  step.  The 
profiles  arc  distinctly  different  from  those  for  the  undisturbed  boundary  layer  shown  in  figure  5. 
The  breakdown  of  the  tke  between  small  and  large  scale  turbulence  is  shown  in  figure  10,  and 
this  is  clearly  dominated  by  the  large  scale  structure  tke  as  distinct  from  the  small  scale  tke.  This 
gives  an  important  contrast  to  the  results  for  the  undisturbed  boundary  layer  that  was  dominated 
by  the  small  sealc  tke.  In  figure  1 1  the  speetra  of  the  vertical  velocity  fluctuations  is  given  at 
three  different  heights  in  the  boundary  layer.  In  this  ease  the  spectra  show  significant  level 
differences  at  different  locations  aeross  the  boundary  layer,  and  the  theory  follows  the  trend  in 
the  measurements,  as  well  as  providing  the  speetral  shape.  In  contrast  the  results  for  the 
boundary  layer  36  step  heights  downstream  of  the  step  (figure  12)  show  that  the  vertical  velocity 
speetra  are  similar  at  all  locations,  similar  to  the  undisturbed  boundary  layer. 

Finally  in  figure  13  wc  show'  the  estimates  of  the  surface  pressure  speetra  for  the  undisturbed 
boundary  layer,  and  the  boundary  layer  at  9  and  36  step  heights  downstream  of  the  step.  Given 
the  good  predictions  of  the  vertical  velocity  one  might  expect  the  surfaec  pressure  speetra  to  be 
well  predieted  as  well,  given  the  relationship  specified  by  equation  (30).  However  the  only 
reasonable  prediction  is  for  the  undisturbed  boundary  layer  at  low  frequencies.  However,  the 
analysis  given  here,  and  the  measurements  have  only  included  the  turbulcnee  in  the  outer  part  of 
the  boundary  layer  and  have  not  considered  the  viseous  sublayer  and  the  log  law  region,  w'hieh 
arc  closest  to  the  wall  and  will  dominated  the  terms  in  equation  (30).  It  is  not  surprising  therefore 
that  the  predieted  surfaee  pressure  speetra,  that  arc  dominated  by  near  wall  effeets,  are  not  well 
predieted  by  this  method. 

3.8  Conclusions 

In  this  study  we  have  made  use  of  the  unsteady  part  of  the  RANS  equations  to  model  the  flow  in 
a  boundary  and  downstream  of  a  forward  facing  step.  It  has  been  shown  that  the  flow  can  be 
broken  down  into  large  and  small  scale  turbulence.  The  large  scale  turbulent  structures  are 
eonveeted  at  their  own  fixed  speed,  while  the  small  seale  turbulence  is  uncorrelated  from  point  to 
point  in  the  vertical  direction,  and  eonveeted  at  the  local  flow  speed.  The  energy  balance 
between  these  two  meehanisms  is  dominated  by  the  small  scale  turbulence  for  the  undisturbed 
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boundary  layer,  but  in  the  case  of  the  flow  down  stream  of  a  step  results  are  quite  different.  At  9 
step  heights  downstream  of  the  step  the  large  scale  turbulence  dominates,  indicating  that  the 
disturbances  generated  by  the  step  are  dominating  the  flow. 

While  the  estimates  of  the  vertical  velocity  compare  well  with  the  measurements,  the  surface 
pressure  spectra  at  different  locations  dowmstream  of  the  step  are  not  as  well  estimated.  The 
explanation  for  this  is  that  the  model  has  not  included  the  inner  layers  where  the  turbulent 
structures  are  expected  to  dominate  the  surface  pressure  fluctuations. 
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Figure  1:  The  Greens  functions  for  a  boundary  layer  flow  at  a  Reynolds  number  of  5xl0\ 
frequency  (0=  1  and  wavenumbers  a=- 1 . 285 ,/M. Left  figure  is  the  Greens  function  for  the  OS 
equation,  right  figure  is  the  interaction  Greens  function  defined  by  equation  (12). 


Figure  2:  The  Greens  function  for  Squires  equation  for  a  boundary  layer  flow  at  a  Reynolds 
number  of  5x1 0\  frequency  ft)=l  and  wavenumbers  a=-1 .285,{h\. 
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Figure  3:  The  Greens  functions  for  a  boundary  layer  flow  at  a  Reynolds  number  of  5x10', 
frequency  ftM).  1  and  wavenumbers  a=-O.I285,ft=\  .Left  Figure  is  the  Greens  function  for  the  OS 
equation,  right  figure  is  the  interaction  Greens  function  defined  by  equation  (12). 
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Figure  4:  The  three  velocity  components  for  a  large  scale  structure  at  <»=10,a=- 12.85,  /fclO  red 
is  v,  black  is  u,  blue  is  vv 
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Figure  5:  The  upper  Figure  shows  the  boundary  layer  mean  velocity  profile  measured  in  the  VT 
Stability  wind  tunnel  at  a  Reynolds  number  of  1 .73x  105 .  the  lower  Figure  shows  the  measured 
turbulent  stresses.  Red  is  v2,  green  is  w2  and  black  is  u2. 
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Figure  6:  The  velocity  mode  amplitudes  for  a  boundary  layer  flow  at  a  Reynolds  number  of 
1 .73x1 05. 
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Figure  7:  The  turbulent  kinetic  energy  budget  for  a  boundary  layer  flow  at  a  Reynolds  number  of 
5xl05  showing  the  breakdown  between  the  energy  in  the  large  scale  structures  and  small  scale 
turbulence. 
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Figure  8:  The  normal  velocity  spectra  for  a  boundary  layer  flow  at  a  Reynolds  number  of 
1 .73x10'  at  three  different  location  in  the  boundary  layer.  Lines  are  experimental  measurements 
and  dots  are  the  model  estimates. 
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Figure  9:  The  upper  figure  shows  the  boundary  layer  mean  velocity  profile  measured  in  the  VT 
Stability  wind  tunnel  9  step  heights  downstream  of  a  forward  facing  step,  at  a  Reynolds  number 
of  1 .73xl05  .  The  lower  figure  shows  the  measured  turbulent  stresses.  Red  is  v2,  green  is  vv2  and 
black  is  u2. 
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Figure  10:  The  turbulent  kinetic  energy  budget  for  a  boundary  layer  flow  9  step  heights 
downstream  of  a  forward  facing  step  at  a  Reynolds  number  of  1 .73x1 05  showing  the  breakdown 
between  the  energy  in  the  large  scale  structures  and  small  scale  turbulence. 
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Figure  1 1 :  The  normal  velocity  spectra  for  a  boundary  layer  flow  9  step  heights  downstream  of  a 
forward  facing  step  at  a  Reynolds  number  of  1 .73x1 05  at  three  different  location  in  the  boundary 
layer,  y/S  =0.0672,  0.2353,0.5379.  Lines  are  experimental  measurements  and  dots  are  the  model 
estimates. 
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Figure  12:  The  normal  velocity  spectra  for  a  boundary  layer  flow  36  step  heights  downstream  of 
a  forward  facing  step  at  a  Reynolds  number  of  1 .73x  10*  at  three  different  location  in  the 
boundary  layer,  y/S  =0.0672,  0.2353,0.5379.  Lines  are  experimental  measurements  and  dots  are 
the  model  estimates. 
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Figure  13:  The  surface  pressure  spectra  for  a  boundary  layer  flow.  Black  line  (measurements) 
and  (predictions)  dots  are  for  the  undisturbed  boundary  layer,  blue  line  (measurements)  and 
(predictions)  dots  are  at  9  step  heights  downstream  of  a  forward  facing  step,  and  red  line 
(measurements)  and  (predictions)  dots  are  at  36  step  heights  downstream.  All  cases  are  at  a 
Reynolds  number  of  1 ,73xl05  at  three  different  location  in  the  boundary  layer,  \/8  =0.0672, 
0.2353,0.5379. 
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